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We present deep wide- field i^^-band observations of six Hyperluminous Infrared Galax- 
ies (HLIRGs) spanning a redshift range 0.44 < z < 1.55. The sample resides in a wide 
variety of environments, from the field to Abell 2 clusters, with a mean galaxy-HLIRG 
clustering amplitude of {Bgu) = 190 ± 45Mpc^ '^^. The range in environments, and the 
\ mean clustering level, are both greater than those seen in local IR-luminous galaxies, 

from which we infer that the range of galaxy evolution processes driving IR-luminous 
galaxy evolution at z > 0.5 is greater than locally, and includes mergers between gas- 
rich spiral galaxies in the field, but also includes encounters in clusters and hierarchical 
^ ■ buildup. The similarity in the range of environments and mean clustering amplitude 

\-{ I between our sample and QSOs over a similar redshift range is consistent with the in- 

C/3 . terpretation where evolutionary connections between IR-luminous galaxies and QSOs 

■ are stronger at z > 0.5 than locally, and that, at these redshifts, the processes that 

drive QSO evolution are similar to those that drive IR-luminous galaxy evolution. 
From comparison of the HLIRG and QSO host galaxies we further postulate that a 
/\ ' larger fraction of IR-luminous galaxies pass through an optical QSO stage at z > 0.5 

^ \ than locally. 

Key v^rords: clusters: galaxies - galaxies: evolution - galaxies: starburst - galaxies: 
active - infrared: galaxies 



1 INTRODUCTION 

Since the discovery by the Infrared Astronomical Satellite 
(IRAS) in 1983 of a large population of galaxies with sig- 
nificant infrared (IR) emission, substantial effort has been 
expended on understanding the nature of the most luminous 
end of this IR galaxy population. These sources, termed 
Ultraluminous Infrared Galaxies (ULIRGs) if their IR lu- 
minosity exceeds 10^^1/0, and HLIRGs if their IR lumi- 
nosity exceeds IQ^^Lq, are found over a very wide range 
in redshift, with mo st lying at z < 0.1 dSoifer et aJll984 : 
ISaunders et a but with a significant number lying in 

the range 0.1 < z < 0.4, and a few lying at higher red- 
shifts. Although the consensus is now that ULIRGs and 
HLIRGs are powered by some combination of violent star 
formation and black hole accretion surrounded by large 

* E-mail: duncan@ipac.caltech.edu 



masses of gas and dust, the triggers for this activity, and 
how these galaxies evolve, are not known. Locally, ULIRGs 
are thought to be mergers between two or more gas-rich 
spiral galaxies, taking place almost exclusively in poor en- 
vironments, and that a small number of these ULIRGs 
evolve into optically selected QSOs dSoife r et al| | l984l: 
Leech et aJll994l:ISanders fc Mirabeillig gG: Rigopoulo u et all 
^99; Fafrah_et_aJl20M J et ali.200Z : .Tacconi et all 

20021: iFarrah et allboOSt) . At higher redshifts however, the 
picture is less clear. It has been suggested that a greater va- 
riety of galaxy formation processes ma y play a role amongs t 
ULIRGs and HLIRGs at high redshift (iFarrah et alll2002bD . 
a change which may manifest itself in their environments. 

The Hyperluminous Infrared Galaxies (HLIRGs) , which 
generally lie ai z ^ 0.3, have been studied extensively 
since their discovery, motivated by their extreme lumi- 
nosities which make them amongst the most luminous 
objects in the Universe. The first HLIRG to be found 
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Table 1. Hyperluminous Infrared Galaxy sample 



Name 


z 


RA 


(J2000) 




Dec 


Spectrum 






Exp. time (s) 






IRAS F00235+1024 


0.58 


00 


26 


06.7 


10 


41 


27.6 


nl 


17.05 


13.15 


1521 


20.68 


16.51 


IRAS P09104+4109 


0.44 


09 


13 


45.4 


40 


56 


28.0 


Sy2 


15.00 


13.24 


2341 


20.81 


15.90 


IRAS F10026+4949 


1.12 


10 


05 


52.5 


49 


34 


47.8 


Syl 


16.85 


14.00 


3511 


21.02 


17.95 


IRAS F10119+1429 


1.55 


10 


14 


37.8 


14 


15 


59.7 


QSO 


16.37 


~ 14.31 


3862 


20.87 


18.69 


LBQS 1220+0939 


0.68 


12 


23 


17.9 


09 


23 


07.3 


QSO 


17.54 


~ 13.08 


1170 


20.53 


16.97 


IRAS F14218+3845 


1.21 


14 


23 


55.5 


38 


31 


51.3 


QSO 


17.21 


13.26 


3277 


21.15 


18.21 



Magnitudes are taken fro m the data presen t ed in this paper. Infrared (1 — 1000/im) luminosities, given in units of bolometric solar 
luminosities, are taken from lRowan- RobinsonI <2000l) and lFarrah et aj ^2002h^ an d reseated to A = 0.7 , = 0.3 and Hq = 70. kg is 
the faintest object detected by SExtractor for each field, mj^ was derived from iPozzetti et aj l2003l) for the redshift of each object. 



(P09104+4109, At z = 0.44) was a cD galaxy in the core 
of a rich cluster, identifie d to have extreme IR emission 
by iKleinmann et aj 1^8^, with a far infrared luminosity 
of 1.5 X IQ^^h'^Q Lq. Then, in 1991, Rowan-Robinson et al 
Jl99ll) identified F10214+4724 at ^ = 2.286, with an appar- 



Later obser- 



ent far infrared luminosity of 3 x 10 /iJq Lq 
vations revealed a large mass of molecular gas (10^^ 
j jBrown fc vanden Bout"l99l': 'Solomon. Downcs fc Radf ordI 
ll^9^V r"a Sevfert emission spectrum ( Elston ct al. 1994), and 
evidence for lensing with a magnification of about 10 in the 
infrared jGraham fc Liulll995l: iB roadhurst fc Lehaijll995t 
lEisenhar^^^^ Il996l : iGreen fc Rowan- Rob inaonl Il99d) . 
These objects appeared to presage a new class of infrared 
galaxy. 

Later observations of larger samples of HLIRGs uncov- 
ered a more detailed picture. Hubble Space Telescope (HST) 
imaging (Farrah ct al 2002a) revealed that a wide range of 
morphologies are present in the HLIRG population, from 
merging systems to QSOs in apparently relaxed systems. X- 
ray, IR, and sub-millimetre observations showed that, in all 
cases, HLIRGs are powered by a mixture of dust-enshrouded 
black hole accretion and violent star formation, with inferred 



star f ormation rates of ^ SOOMpyr ^ jRowMbRobmsoE 
20001: IVerma et a] l2002l: iFarrah et al l2002bl: IWilman et a 



20031) , suggesting that HLIRGs are comprised of both merg- 



ers between gas-rich spiral galaxies, and young galaxies go- 
ing through their maximal star formation epochs whilst har- 
bouring an AGN. 

Despite this progress, the role of HLIRGs in the broader 
picture of galaxy and AGN evolution remains unclear. It 
is not known whether HLIRGs as a class are a simple ex- 
trapolation of the local Ultraluminous Infrared Galaxies 
(ULIRGs, Lir > 10^^ 1/0 ) making them mostly mergers be- 
tween gas-rich spirals, or whether a wider range of galaxy 
formation processes play a role in HLIRG evolution. Also, 
the links between HLIRGs and QSOs at comparable red- 
shifts are not well understood. Locally, it is thought that 
some fraction of ULIRGs evolve into optically selected QSO s 
JSanders et aJll988l: iFarrah et a'ill200ll : iTacconi et a ? 2002) , 
but it is not known whether this is also true in the distant 
Universe. 

Many of these unknowns result from two major ob- 
stacles in studying HLIRG evolution. Firstly, HLIRGs con- 
tain very large masses of gas and dust, making observations 
of the galaxies themselves at all wavelengths (except per- 
haps the far-infrared and sub-millimetre) prone to obscu- 
ration bias. Secondly, the presence of a luminous starburst 



and AGN in all HLIRGs means that observations will be 
affected by the orientation of the HLIRG relative to us. 
These problems can however be partly overcome by exam- 
ining the environments of HLIRGs. Since the determina- 
tion of environments is independent of orientation and dust 
content, they are a useful tool in studying AGN evolution, 
and have been used extensively in studying both normal 
and active galaxies (Lonaair fc Scldnor 1979; Yoo fc Green] 
1987 : Hill fc Lillv 1991: Loveday et al.l995; Wold ct al 2P.0a 
2001 ilMcLure fc Dunlori200 jj^lSanchez fc Gonzalez-Serranol 



2002). Studying the environments of HLIRGs therefore can 
help clarify the relations between HLIRGs and other AGN 
classes. 

In this paper, we investigate the environments of six 
HLIRGs, using deep wide field Ks-hand imaging. Observa- 
tions are described in §2 and analysis is described in §3. Re- 
sults are presented in §4, with discussion in §5. Finally, our 
conclusions are summarized in §6. Unless otherwise stated, 
we assume A = 0.7, Qo = 0.3 and Ho = 70 km s"^ Mpc"^ 



2 OBSERVATIONS 

We selected for observation six HLIRGs from the sample pre- 
sented by Rowan- Robinson (2000). The sample, their red- 
shifts and other basic data are presented in Table Five 
of these objects were selected to lie approximately in the 
redshift range 0.6 < z < 1.6, where the greatest evolu- 
tion in the IR galaxy population is thought to occur (e.g. 
iRowan-Robinson et a] lll997l) l. Additionally, we observed 
one further HLIRG, P09104-I-4109, which lies a.t z = 0.44 
and is already known to lie in a rich cluster, to act as a con- 
trol for our observation and analysis methods, though we do 
not include this source in the discussion. None of our targets 
show any evidence for significant gravitational lensing. 

Observations were made on 25-26th December 2001 us- 
ing the INGRID wide field near-infrared imager and a Ks- 
band filter, on the 4.2m William Hershel Telescope (WHT). 
INGRID is a 1024 x 1024 pixel array, with a scale of 0.238" 
pix~^, corresponding to a field-of-view of ~ 17' . At the red- 
shifts of our sample this corresponds to a physical field-of- 
view of > 1.36Mpc. The targets were centred approximately 
in the INGRID field-of-view, with exposure times selected 
to reach a minimum depth of M^^ -I- 2 at the redshift of each 
object. Observing conditions were generally good, with little 
cloud cover and seeing of ~ 0.8", however the atmospheric 
stability was variable, resulting in non-photometric nights. 
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The total exposure time for each object was divided into 
several nine point 'box' dither patterns, with a 16.7" offset 
between each position, to allow the subtraction of cosmic 
rays, hot pixels, and the infrared sky background. As the 
INGRID field of view is too small to reliably estimate field 
galaxy counts from the edges of the HLIRG fields, separate 
control fields were also observed, with similar galactic lat- 
itudes and exposure times as the sample. For photometric 
calibration we observed a selection of infrared standard stars 
throughout each night, at several different airmasses. 



3 DATA REDUCTION 

Following debiasing and fiatfielding, the data were reduced 
using our own custom-written IRAF pipeline, based in part 
on the Quicklook INGRID data reduction pipeline from the 
Isaac Newton Group (ING). As the exposure time at each 
position in the dither pattern was only 30 seconds, partic- 
ular care was taken in accounting for bright sources when 
subtracting sky noise. Source masks were created for the 
individual frames in each dither pattern by first creating 
an approximate estimate of the sky by median combining 
the nine frames in each dither pattern without applying the 
dither offsets. This 'dummy' sky frame was then subtracted 
from each frame in the dither pattern to reveal the brightest 
sources, which were then masked out when creating the 'real' 
sky frame. The nine frames in each dither pattern were com- 
bined using the IRAF task IMCOMBINE, with dither offsets 
calculated by centroiding two or more bright sources com- 
mon to each frame. After each dither pattern was combined 
into a single image, sky subtraction was performed using 
the 'real' sky frames and pixel masks described earlier. It 
was found that this sky subtraction method worked well in 
dealing with variations in the infrared sky between different 
dither patterns, over the long exposure time for each object. 
The combined images from each dither pattern were then 
stacked together to produce a final image /: 



(1) 



where U and ai are the individual images from each dither 
pattern, and their standard deviations, respectively. The 
standard deviations were derived from the noise fiuctuations 
in each image. Noisy edges, where the total exposure times 
were shorter due to the dither pattern, were clipped off. 
The resulting field of view, at 13' x 13', was still sufficient to 
quantify the environments of the sample. The final images 
were of excellent quality, and fiat to better than 1% across 
the width of the frame, and were in all cases of much higher 
quality than the images from the Quicklook reduction. Our 
reduction pipeline, 'INREP', is available for general use via 
the Isaac Newton Group (ING) web pages ^ . 

Sources we re extracted and catalog ued using the SEx- 
tractor package iBertin fc Arnoutsll996l) . For source extrac- 
tion, we adopted the conservative criterion that a source con- 
stitutes at least four contiguous pixels, with a significance of 
detection of at least 3cr above the background. The default 
SExtractor extraction filter was used to detect faint extended 



http: / /ing. iac.es/Astronomy/Ingrid 



objects, with 32 de-blending thresholds, a cleaning efficiency 
of 1 and a contrast parameter of 0.005. The background es- 
timation was mesh based, with a mesh size of 64 and a filter 
size of 3. As many of the sources in the frames were faint 
and slightly extended, the magnitudes were calculated using 
corrected isophotal (ISOCOR) magnitudes within SExtrac- 
tor. The calibration zeropoint magnitudes for each object 
frame and control frame were calculated using the observa- 
tions of the infrared standard stars. The errors on the final 
magnitudes arising from the standard star calibrations are 
Am = 0.05, with a further error of Am = 0.05 from the 
source extraction. To correct our measured magnitudes to 
the rest-frame J^s-band, we computed fc-corrections assum- 
ing a standard power law, with a = 2. 



4 ANALYSIS 



To evaluate the environmental richness of each object in 
our sample we use two independent clust ering measures; the 
-BgE galaxy-galaxy correlation statistic jLongair fc Seldneil 
which for our sample we refer to as the Bgh 
:alaxy-HLIRG cor relation statistic, and the No.^ statistic 
HiU fc Lilly! Il99j) . The Bgu statistic is the amplitude of 
the spatial cross-correlation function, and relies on knowing 
the volume density of galaxies around the HLIRG, which it- 
self requires knowledge of the luminosity function (LF) . The 
Nq.z statistic is a simple, more direct counting statistic, and 
involves counting all the sources within 0.5Mpc of the target 
that lie within a certain magnitude range. Although the Bgh 
statistic is physically more meaningful than the statis- 
tic, the near-infrared luminosity function is not known with 
great precision at high redshifts. Therefore, we have adopted 
the Bgh statistic as our primary measure of environmental 
richness, but we also use the statistic as an independent 
check. 

To compute the _Bgh sta tistic, we follow the prescription 
of lLongair fc Seldneil ^ll979^ . As this statistic is not trivial to 
evaluate, we provide a detailed outline of the procedure here. 
The spatial cross-correlation function, £,(r), can be described 
as a simple power law: 



and is derived from the 
correlation function w(9): 

w{e) = A^i,e'-^ 



(2) 

galaxy-HLIRG angular cross- 



(3) 



by deprojecting from the celestial sphere. The amplitude of 
the angular cross- correlation function, Agh, is essentially the 
richness measure of the number of galaxies about a point on 
the sky, and is converted to -Bgh, the richness measure of the 
number of galaxies about a point in space, by translating the 
angular cross-correlation function w{9) to the spatial cross- 
correlation function f (r). 

To perform this conversion, the Agh statistic is evalu- 
ated: 



Agh = 



Nt ~ Nb 3 - -y 
Nb ~ 



(4) 



where Nt is the total number of galaxies within a radius 
8 around the target, corresponding to 0.5 Mpc at the tar- 
get redshift, and Nb is the number of galaxies within the 
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Table 2. HLIRG clustering statistics 



Name 














Abell class 








xlO^rad"'^^ 




Mpci-''^ 






IRAS F00235+1024 


7.69 


2.77 


0.66 


0.28 


208.63 


88.31 


field 


IRAS P09104+4109 


29.65 


5.45 


2.75 


0.40 


819.75 


119.77 


2 


IRAS F10026+4949 


17.81 


4.22 


1.73 


0.40 


611.46 


141.01 


1 


IRAS F10119+1429 


8.88 


2.98 


0.72 


0.31 


384.32 


167.33 





LBQS 1220+0939 


3.77 


1.94 


0.40 


0.26 


143.78 


94.25 


field 


IRAS F14218+3845 


4.45 


2.11 


0.06 


0.21 


21.84 


81.27 


field 



Quoted errors only include the counting errors. The methods used to calculate the Bgj, and A^o.s statistics, and the Abell classes, are 

described in Jl]and JF] respectively. 



same radius in t he control field for tliat object. We make 
the assumption llWold et al 2000, 200 J) that 7 = 1.77. 
The precise choice of 7 will however not affect the de- 
rived clustering amplitudes , as it has previously been shown 
ifPresta ge fc PeacocklllQS^ . Il989^ that Bgh is insensitive to 
the choice of 7, as long as 7 ~ 2. 

The Bgh statistic can now be calculated, and normalised 
to the integral LF, $(m;, z) which represents the number of 
galaxies more luminous than mi per unit co-moving volume 
at redshift z: 



^{mi,z)L, " 



(5) 



where pg is the average surface density of background galax- 
ies, de is the angular diameter distance to the target and J-y 
is an integration constant: 



r^[(7 + i)/2] 
r(7) 



3.78 



(6) 



The integration to <^( mi,z) is perform ed by taking a uni- 
versal LF of the form (ISchechterlll97el) : 



c^{L) ~ (07L*) [L/L'T exp {-L/V) 



(7) 



and integrating down to the completeness limit at the target 
redshift L{mi, z): 



Mmi,z) 



(P{L)dL 



(8) 



L(mi,z) 



For the parameters used in the LF, we use the most recent 
determination of th e high redshift near-IR LFs as given by 
iPozzetti et all ll2003l) . who have determined the evolution of 
the near-IR LF in the J- and Jfs -bands in redshift bins of 
z — (0.20,0.65) and z — (0.75, 1.30), using a spectroscopic 
survey of a magnitude limited sample of galaxies with Ks < 
20. In terms of absolute magnitudes the luminosity function 
can be expressed as: 



0(Af) = C/(M)('^+^) exp(-/(M)) 

where the constant C = 0.4 ln(10)(j!)* , and: 

/(A/) = io(-«'^(^^-^'^*» 



(9) 



(10) 



As some of the galaxies i n our sample li e outside of the 
redshift bins given by Poz zetti et all (|200^, we must make 
assumptions as to the nature of the LF at these magnitudes. 
In the redshift range 0.2 < z < 0.65, we have assumed a = 
-1.25, M* = -25.64 and <;/>* = 6.11 x 10"*Mpc"^ In the 
redshift range 0.75 < z < 1.55 we have assumed a = —0.98, 



M* = -2 5.54 and <?!)' = 9.98 x lO-'Mpc'^ We note that 
IPozzettret al ( 2003) find only mild evolution in the J^s-band 
luminosity function over the redshift range < z < 1.3, 
hence our use of their luminosity function at z — 1.55 is 
unlikely to be a major source of error. 

The errors in Agh and B^i^ were calculated following the 
prescription given bv lYee fc Lopez-Cruj lll999f) : 



A 



gh 



B. 



gh 



y^{NT-NB) + 1.3^NE 

Nt - Nb 



(11) 



which is based only upon (non-Poissonian) counting statis- 
tics. Other systematic errors are discussed in H5.2I 

To com pute the statistic, we follow the procedure 
described bv lHill fc Lilhl il99 J) . This statistic is computed 
by counting all sources within a 0.5 Mpc radius of the tar- 
get, with a magnitude in the interval (m, m + 3), where m 
is the magnitude of the target. Subtracted from this number 
is the expected number of background field galaxies in the 
same magnitude interval calculated from the control fields. 
The resulting number is the N0.5 statistic. This measure- 
ment does not require integration of a luminosity function, 
although if the luminosities of HLIRGs and cluster galaxies 
evolve very differently with redshift then the A'o.s statistic 
will give inaccurate results. There is no evidence why this 
should be a major effect however, so we do not take it into ac- 
count. From Table the limiting magnitude for each object 
is fainter than tuhlirg + 3 except for two cases; the image 
for F10119+1429 reaches ruHLiRG + 2.18, and the image for 
F14218-I-3845 reaches mnLiRG + 2.94. In both these cases 
however the limiting depth is sufficiently near mHLiRO + 3 
that the statistic is still usable. 



5 RESULTS 

5.1 Clustering statistics 

Images of the fields around each HLIRG are presented in 
Figure^ In Figure|5|we plot the galaxy counts in a 0.5Mpc 
region around each HLIRG, and the counts for the asso- 
ciated control field. The Bgh and A'o.s statistics for each 
object, and their errors, are given in Table |5| We also quote 
the approximate Abell classes of our sample in this table, 
based on the conversions given by iHill fc Lilly ( 199lll , in 
our cosmology. We note however that these conversions are 
arbitrary, and hence we base our quoted Abell classes on 
both the Bgh and iVo.5 statistics. 

A wide variety of environments can be seen amongst 
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Figure 1. i^s-band images of tlic fields around each HLfRG. An arrow indicates the H^LIRG. Large tick marks correspond to 25". 
L-R: (top row) F00235+1024 and P09104+4109, (middle row) F10026+4949 and F10119+1429, (bottom row) LBQS1220+0939 and 
F14218+3845. 
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Figure 2. Galaxy counts in a circular region 0.5Mpc in radius around each HLIRG. The solid line shows the counts in the field around 

the HLIRG, and the dashed line shows the counts in the control field. L-R: (top row) F00235+1024 and P09104+4109, (middle row) 
F10026+4949 and F10119+1429, (bottom row) LBQS1220+0939 and F14218+3845. Error bars have been omitted for clarity. 
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the objects in the sample. P09104+4109 resi des in an Abell 
2 clu ster, in agreement with previous results jHines fc Will j 
Il99,^ ^. For the remaining objects, three reside in poor en- 
vironments, and two, F10026-I-4949 and F10119+1429, re- 
side in clusters, of Abell class 1 and respectively, though 
we note that the detection of clustering for F10119-I-1429 
is marginal, at 2.5a. Excluding P09104-I-4109, the error 
weighted mean clustering amplitude for the remaining 5 ob- 
jects is (Bgh) = 190 ±45Mpc^ ". 

We can compare the values of and Bgh to each 

other, as these quantities have been used by many previous 
autho rs, and a wel l defin ed relation has been found between 
them. iHill fc Lillvl Jl99 J) derive Bgh oc SOiVo.s. In Figure 01 
we plot this relation, together with Bgh vs. for the 

objects in our sample, a nd the conversion s between Bgh 
and Abell class given bv iHill fc Lillvl lll99lll . We also plot 
the best fit linear relation between Bgh and for our 

da ta. The best fit is well matched to the relation derived 
bv lHiU fc Liiwl lll99ir) and the objects in the sample follow 
this relation closely, hence we conclude that our computed 
values of A'0.5 and Bgh are reasonable. 

5.2 Error budget 

The errors quoted in Table|21for the Bgh and A'0.5 statistics 
are only the counting errors, and do not include three fur- 
ther, potentially important sources of error. In this section, 
we discuss these three error sources in turn. 

The first of these is the luminosity function assumed 
in calculating the Bgh statistic. Due to the uncertainties in 
the high redshift Ka-hand LF used in this paper, this is an 
important issue to address. lYee fc Lopez- Cru3 1^9^ have 
examined the sensitivity of -Bgh to the form of the LF, and 
find that an error in the assumed value of M* of up to 
±0™.5 will still yield essentially the same results. Similarly, 
an error in a of up to ±0.3 will only affect Bgh by at most 
~ 20%. Although the K band LF is not well constrained at 
the redshifts of our targets, we do not expect that the true 
LF differs from our assumed LF by such gross margins, and 
therefore we conclude that uncertainties in our computed 
Bgh statistics due to the assumed LF are at most 10%. 

The second of these error sources is the method used to 
compare the data to the luminosity function in calculating 
the Bgh statistic. This can be done in two ways. The data 
can be fc-corrected to the rest-frame Ks-hand, and compared 
to the Ks-hand luminosity function given bv lPozzetti et all 
l|20o3), or the data can be compared directly to the J band 
luminosity function presented by Pozzctti ot al (2003), with- 
out applying a fc-correction, as observed frame ifs-band ap- 
proximately samples rest-frame J band at the redshifts of 
our sample. Both methods have advantages and disadvan- 
tages; the rest-frame J^s-band is a less contaminated tracer 
of evolved stellar mass than the J band, but suffers from the 
extra uncertainties introduced by applying fc-corrections. Al- 
though we chose to fc-correct our data to the rest-frame Ks- 
band, we also examined the effect of comparing our data di- 
rectly to the rest-frame J-band luminosity function, without 
applying fc-corrections. We computed Bgh and crs ^ with- 
out applying fc-c orrections, and using the J-band luminos- 
ity function from iPozzetti et"ail (|2003), and found that the 
Bgh values were within la of the values computed using 
fc-corrections and the A's-band luminosity function for all 
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Figure 3. Bgh vs. A'0.5 for o ur sample. F i lled he xagons are the 
Abell calibrations taken from iHill fc Lillvl il99ll) . stars are our 
data. The solid lin e shows the relation Bgh oc SOA'o.S found by 
iHill fc Lillvl <199ll) . The dashed line is a least squares fit of our 
data, which gives Bgh = 29.07A'o.5 + 14.95 



the objects except F10119-I-1429, which was within 1.5a. 
Furthermore, each object still resided within the same type 
of environment as before, with P09104-I-4109, F10026-I-4949 
and F10119+1429 residing in clusters (with the same signifi- 
cance of detection of clustering) and the other 3 objects lying 
in the field. The values of Bgh and as ^ were however in all 
cases slightly higher using the J-band luminosity function 
and no fc-corrections, with F10026-I-4949 and F10119-I-1429 
predicted to (just) lie in Abell 2 clusters. We conclude that 
using the Ks-hand luminosity function with fc-corrections, 
or just the J-band luminosity function, will not significantly 
change our results, and that our choice of using the Jfs- 
band luminosity function with fc-corrections was the most 
conservative. 

The third of these error sources is the choice of cosmol- 
ogy. Whilst both Bgh and A'^0.5 are relatively insensitive to 
the choice of A and Qq, the sensitivity to Ho is marked, par- 
ticularly in the range 65 < Ho < 80. This is illustrated in 
Figure 21 where we have plotted Bgh as a function of Ho for 
F10026+4949. By varying Ho over a relatively small range, 
the effect on Bgh is dramatic, changing from Bgh ~ 350 for 
Ho ~ 65 to Bgh ~ 800 for Ho ^ 75. Most current measure- 
ments of Ho produce values in the range 70 to 75, albeit 
with a significant error (Frcodman ct al (2001) and refer- 
ences therein). Our adopted value of Ho = 70, and therefore 
our derived clustering amplitudes, are therefore conserva- 
tive. 

In summary, we conclude that, whilst there are fur- 
ther, significant sources of error than the Poisson errors 
on the derived clustering statistics, none of these sources of 
error should significantly change our results. Furthermore, 
we have in all cases adopted the most conservative method 
possible in computing the clustering amplitudes, making it 
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highly unlikely that any of our objects that we infer to lie 
in clusters actually reside in the field. 



6 DISCUSSION 

6.1 The hyperluminous phenomenon at low and 
high redshift 

Our sample of five objects (excluding P09104), although 
small, exhibit a wide variety of environments, from poor 
environments to Abell 1 clusters. This strongly suggests 
that there is a wide variety of environments amongst ob- 
jects with high levels of IR emission at z ~ 1 generally, 
although a larger sample would be required to confirm 
this. This variety of environments is not seen amongst local 
ULIRGs, which are genera lly not found in rich environments 
iSanders fc MirabeJliooil . The mean environmental rich- 
ness of our sample, at (-Bgh) ~ 190±45Mpc^'^^, is higher, at 
just over 3a significa nce, than the galaxy-galaxy correlation 
statistic both locally iLovedav et ail 995: Guzzo ct al 1997D 
and at moderate redshifts H^^doi^c Lilly... 1996 ') , though we 
note that the latter study measured values of Bgg up to red- 
shifts of only z ~ 0.5. We infer that, in going from z; ~ 
to ~ 1 the environments of IR-luminous galaxies become 
more diverse, and that the mean environment becomes richer 
than that of normal galaxies both locally and at moderate 
redshifts. This is supportive of the idea that a wider va- 
riety of galaxy formation processes are important amongst 
the IR-luminous galaxy population at high redshift than lo- 
cally, such as hierarchical buildup or encounters in clusters. 
Indeed, one object in our sample, IRAS F10026-I-4949, may 
be an example of such a galaxy. This sourc e harbours both a 
starburst and an AGN (iFarrah et ai20 02b'). and lies in a rich 
cluster. From HST imaging, this source also pos sesses mul- 
tiple very close companions jFarrah et a^l2002a^ . It is thus 
an excellent candidate for being a cD galaxy in the process 
of formation in a cluster at z = 1.12, and may be the higher 
redshift analogue of the clustered IR-l uminous active galax- 
ies P09104-h4109 and P18216-h6419 llSchneider et allll992l: 
iHines fc Wiii3ll993tlWold et aJl2002l) . 



6.2 Comparison to Quasar environments 

Since environments are independent of orientation, we can 
compare the environments of our sample to those of other 
classes of active galaxies at comparable redshifts to ex- 
amine possible relationships between HLIRGs and other 
AGN classes. The environments of AGN over a wide red- 
shift range have been studied by many authors in efforts to 
disentangle the myriad AGN taxonomy. Radio Loud QSOs 
(RLQs) and Radio Galaxies (RGs) are found in a diverse 
range of environments, from the field to Abell Class 2 and 
greater, at both low and high redshift ([L ongair fc Seldneil 
19791: lYee fc Greenlll987t iPrestaee fc Peacock! Il988l | l989|: 
Hill fc Li llv"l99l': 'AUington-Smith et al"l993': 'Zirbei ll997l: 
Wold et a l 2000; Sanc hez_fc G onzalez-Serrano 200^. Over- 
all, RLQs and RGs appear to prefer moderately rich envi- 
ronments on average, of around Abell Class 0. Determining 
whether or not there is evolution in the environments of 
RLQs or RGs with redshift is difficult, due to gravitational 
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Figure 4. The clustering amplitude of F10026-I-4949 as a func- 
tion of Ho (for A = 0.7, Co = 0.3). Error bars are the Poisson 



lensing and selection biases, but from recent results it ap- 
pears that there is no significant difference in RLQ and RG 
environments at 2 < 1.0 (Wold ot al 2000). For Radio Quiet 
Quasars (RQQs) a broadly comparable picture has now 
emerged. RQQs are found in a similarly diverse range of en- 
vironments to RLQs, from the field to Abell class 2 or richer, 
over a wide redshift range fYee & Green '1987"; DunloD ot all 
[l993; Fisher ctal ,1996 ; Dehorn ot al 1997: .Tanak a ot a| 
^OOOt iMcLure fc Dunloa 2001). Indeed, the most recent re- 
sults, based on deep optical imaging, show that the envi- 
ronments of RQQs and RLQ s at moderate r edshifts are 
statistically indistinguishable jWold et "5I I2OOJ) . At z ^ 1 
there is as yet no clear consensus, but the environments of 
RLQs and RQQs at these redshifts do not appear to be 
significant ly different ij Hintzen^ Rpmanishin fc Valde slll99ll : 
Fllutchings . Crampton fc^ohngonl 1995: Hutchings 19951). 

In the local Universe, it was initially suggested that 
ULIRGs as a class are precursors to optically selected QSOs, 
although later results show that this is pro bably true for 
only a small subset of the ULIRG population fe anders et all 
•l988; Farrah ct al 2001; Taccorfi ct al 2002). This is reflected 
in their environments, as locally ULIRGs and QSOs reside 
in different environments, with ULIRGs generally lying in 
the fleld and QSOs lying in moderately rich environments 
on average, with a diverse range. We can compare the range 
of environments in our sample of HLIRGs to those of QSOs 
at comparable redshifts to see if this is also true at high red- 
shift. The diverse range of environments seen in our sample 
of HLIRGs is qualitatively similar to the range seen in RLQs 
and RQQs, and our mean value for Bgh is comparable to the 
mean galaxy-quasar correla tion statistic. Bag, for R LQs and 
RQQs at similar redshifts J Wold et aJl200Ct |2001^ . This is 
illustrated in Figure |S] where we plot clustering amplitudes 
vs. redshifts for our HLIRGs, together with a representative 
sample of clustering amplitudes for QSOs from the litera- 
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ture iYee k, Ellingsonlll99ilHall fc Greenlll998l: IWold et all 
l2000l . l200iri . converted to our cosmology. The error bars on 
the QSOs have been intentionally omitted for clarity, but 
are comparable in size to the error bars on the HLIRGs. 
We note that the clustering measures were derived using 
different data, with our clustering amplitudes derived from 
Ks-hanA data and the clustering amplitudes from the other 
samples derived from multiband optical data, however it is 
unlikely that this will introduce any systematic biases be- 
tween the datasets. Therefore, we extrapolate from this that 
evolutionary links between IR-luminous galaxies and opti- 
cally selected QSOs may become stronger with increasing 
redshift. Furthermore, we postulate that, at approximately 
z > 0.5, the range of galaxy evolution processes driving 
IR-luminous galaxy evolution caused by the greater diver- 
sity in environments is similar to the range driving QSO 
evolution. We note that one object in our sample, LBQS 
1220-1-0939, is a QSO from the Large Bright Quasar Survey 
jHewett. Fohz fc Chaffe3ll995l) that was only later discov- 
ered to be hyperluminous in the IR; the selection criteria 
for this object therefore make it unsuitable for determin- 
ing evolutionary connections between HLIRGs and QSOs. 
Even with LBQS 1220-f 0939 removed however, the sample 
still contains the same diverse range of environments, and 
the mean clustering amplitude is still comparable to that of 
QSOs at similar redshifts. 

We can take this one step further via a comparison 
of HLIRG and optically selected QSO host galaxies. QSOs 
at all redshifts are thought to lie in luminous, > L* host 
galaxies, with elliptical hosts becoming more prevalent for 
both R LQs and RQ Qs with increasing optical nucl ear lumi- 
nosit y iMcLure et al 1999 ; Percival ot al 2001 : Dunlop et all 
|20o3). Conversely, although local ULIRGs are thought to 
evolve into ellipticals with a few also passing though a QSO 
phase, are thought to be mergers between (optically) sub- 
L* galaxies, w ith the majority o f ULIRGs themselves also 
being sub-L* dColina et aJ 120011) . For the HLIRGs in our 
sample, this is not the case. One object, F14218-I-3845, is 
a QSO with a host magnitude that suggests (similarly to 
the other QSOs in t he sample) the hos t is a very lumi- 
nous, massive system jFarrah et a^l2002a^ . For the two non- 
QSOs in our sample, the picture is less clear, with one object 
(F00235-I-1024) being 0.5 magnitudes dimmer than mj^-, but 
the other (F 10026-1-4949) over a magnitude brighter than 
m^. We note however that the J-band magnitudes of those 
HLIRGs without a QSO presented bv lFarrah et aJ ll2002afl 
are, with the exception of F00235-I-1024, much brighter than 
m*i. Whilst more data is evidently needed, it appears rea- 
sonable to assume that IR-luminous galaxies may become 
more similar to QSO hosts with increasing redshift. Whilst 
this reinforces the idea that evolutionary links between IR- 
luminous galaxies and QSOs become stronger with increas- 
ing lookback time, it also suggests that, at 2 ^ 0.5, a greater 
number of IR-luminous galaxies evolve directly into optically 
selected QSOs than locally. 



T CONCLUSIONS 

We have presented deep wide-field K-band imaging of the 
fields of six Hyperluminous Infared Galaxies, and quantified 
their environments using the Bgh galaxy-HLIRG correlation 
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Figure 5. The clustering amplitude vs. redshift for the 
HLIRGs in our sa mple. Also plotted are cluster ing amplitudes fo r 
QSOs taken fromlYee fc Ellingsonl ir993'). Ha ll fc GreenI <1998l) . 
IWold et all i20n(t and IWold et alnOOiS) . converted to our cos- 
mology. The HLIRGs are plotted as larger symbols for clarity. 



amplitude, and the clustering statistic. We conclude the 
following: 

1) The HLIRGs in our sample reside in a diverse range of 
environments, from the field to Abell 2 clusters. The mean 
clustering level of the sample, at {Bgh) ~ 190 ± 45, and the 
range of environments, are both significantly greater than 
those of the most luminous IR galaxies locally. We infer 
that, at high redshift, the galaxy evolution processes driv- 
ing the evolution of IR-luminous galaxies are more diverse 
than at low redshift, and include mergers between gas-rich 
spirals in the field, but also include encounters in clusters 
and hierarchical buildup. 

2) The mean clustering amplitude of the sample, and the 
range in environments, are comparable to those of QSOs 
over a similar redshift range. We postulate from this that, 
at z 0.5, the range of galaxy evolution processes driving 
IR-luminous galaxy evolution is similar to the range that 
drives QSO evolution. When combined with the similarities 
between HLIRG host galaxies and QSO host galaxies at the 
same redshifts, this further suggests that a greater fraction of 
IR-luminous galaxies evolve directly into optically selected 
QSOs at high redshift than do locally. 
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